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In this work, we present a topological characterization of superconductivity in a prototype electron 
fractionalization model for doped Mott insulators. In this model, spinous and holons are coupled via 
the mutual Chern-Simons gauge fields. We obtain a low-lying effective description of the collective 
current fluctuations by integrating out the matter fields, which replaces the conventional Ginzburg- 
Landau action to describe the generalized rigidity of superconductivity. The superconducting phase 
coherence is essentially characterized by a topological order parameter related to a Gaussian linking 
number, and an experiment is proposed to probe this topological property. We further show that a 
gauge-neutral fermionic mode can naturally emerge in this model, which behaves like a Bogoliubov 
quasiparticle. 

PACS numbers: 74.40.Kb,74.72.-h 

I. INTRODUCTION 

The microscopic nature of high-Tj, superconductivity in the cuprates is still under heavy debate after more than 
two decades of intensive studies. It has been widely accepted that the cuprates may be properly considered as doped 
Mott insulators and the superconductivity arises after the antiferromagnetic (AFM) long range order is destroyed by 
doping. The challenge is that a typical doped Mott insulator is a strongly correlated electron system in which two 
electrons staying at the same lattice site will cost a huge energy as compared to other energy scales like the hopping 
energy. With the failure of conventional perturbative many-body methods, to tackle a doped Mott insulator one 
usually resorts to either the so-called Gutzwiller projected trial wavefunction or constructing an effective low-energy 
theorjff'. In the latter approach, the electron fractionalizatiorP^ is often utilized in order to implement some key 
features of a doped Mott insulator. Indeed, the fractionalization due to the Mottness is most transparent at half- 
filling where only the spin degrees of freedom of the electrons are left, with the charge part being totally frozen out. 
Only by doping can charge carriers emerge which account for a small portion of the total number of the electrons at 
low doping. 

The real challenge is that there exist many different ways of introducing fractionalized elementary objects to 
implement the same local Mottness, i.e., the no double occupancy constraint. These elementary particles are generally 
further coupled to some emergent gauge field(s) due to the uncertainty in the fractionalizatiorP. In the lack of precise 
energetic comparison, a full self-consistent development of a potential fractionalization/gauge theory becomes essential 
in order for it to be convincingly verified or falsified by experiment as well as by general theoretical considerations. 
In principle, an inappropriate fractionalization would always lead to strong gauge fluctuations rendering the effective 
theory intractable. Therefore, a stable, well-controlled gauge theory description will be a highly desired candidate for 
the correct low-energy theory of cuprate superconductors. 

We will consider a particular electron fractionalization description with a mutual Chern-Simons gauge structurd^H^. 
Specifically, an electron in this theory is "fractionalized" into an S — 1/2 bosonic spinon and a charge +e bosonic 
holon, which are minimally coupled to U(l) gauge fields A'^ and A", which do not have their own dynamics, but are 
"entangled" by a mutual Chern-Simons topological term 

Cues = -e'^-^A^^d.A'',. (1) 

which is (2-|-l)D realization of the topological BF theory.'^l Physically it simply means that the spinous and holons 
will perceive each other as quantized tt fluxoidJ^^. In the limit of zero hole concentration, the spinon condensation 
can naturally lead to an AFM long-range ordered state. At finite doping, the mutual Chern-Simons term in Eq. ([T]) 
captured^ a precise "statistical sign structure" identifiecP^I^ in the t-J model. By the tt fluxoids carried by the mobile 
holons, the AFM state will be destroyed and eventually a superconducting will set in with the holon condensation, 
and a zero-temperature phase diagram as a function of doping concentration has been recently discussed in Refs. 
I6I7I Finally, there are other quantum systemi^SHIIl -which have the mutual Chern-Simons gauge structures in different 
contexts. 
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In this paper, we focus on the superconducting state to develop an effective description based on the mutual Chern- 
Simons gauge theory in Sec. II. Instead of a conventional Ginzburg-Landau theory in description of the fluctuations 
of the pairing order parameter, we obtain a gauge description of the low-lying collective fluctuations of charge and 
spin currents, which are essentially the phase fluctuations of the superconducting order parameter. We show that the 
Meissner effect in such a superconductor is basically controlled by a spin gap. In Sec. Ill, we further identify a non- 
BCS-type order parameter I that characterizes the SC phase coherence. Here I is related to the correlation function 
of two Wilson loops (of A" and A'^ ) and is topological as it merely depends on the Gaussian linking number under a 
continuous deformation of the loops. Such a relation is similar to the connection between the Chern-Simons theory 
and Knot theory first studied in Ref. \W[ But by coupling to the matter fields, we show that the spin correlations 
strongly renormalize I such that the phase coherence terminates at a critical point where the spin gap vanishes. In 
Sec. Ill C, we shall propose a direct experimental probe of such a topological nature of the gauge interaction between 
the spin and charge degrees of freedom. Finally, in Sec. Ill D, we show that a fermionic excitation can emerge in the 
superconducting phase, which is gauge neutral and thus coherent, resembling a Bogoliubov quasiparticle. 

II. EFFECTIVE FIELD THEORY DESCRIPTION OF SUPERCONDUCTIVITY 

A. Lattice mutual Chern-Simons theory 

The lattice field theory formulation of the phase string theorj^ for a doped Mott insulator has been recently 
constructed in Refs. I6l7l in terms of partition function: 

Z= [ D[A':A^]D[h\h,b\b]e-' (2) 

in which the action S = C {x denotes spacetime coordinates) with 

C ~ Ch + Cs + Cmcs ■ (3) 

Here the first two terms, Ch and C^, are the Lagrangian densities for the matter fields, holon hi and spinon bi^, 
which carry charge and spin degrees of freedom, respectively. They are given by 

£h = hlido - tA'„ - lA'o + A")/i, -thJ2 {e'^'-^'^^-hlh^^a + h.c.) + ^ [hlh^Y , (4) 

= Y: bl {do - ^aA'i + A^) 6.. - JsJ^ {^^^^bl+aA-. + ^-c) + f (E ^1^.) ■ (5) 

Equations Q and ^ show that the bosonic holon and spinon fields minimally couple to two lattice U(l) gauge fields, 
and respectively, in 2+1 dimensions (/i — a,0, with a — x, y). {A^^ is the external electromagnetic field which 
solely couples to the holons^.) The parameters th and Js are effective hopping integral and superexchange coupling, 
respectively, and do is the difference operator defined along the imaginary time axis. The Lagrangian multipliers, 
A'*'^, control the total numbers of holons and spinous, and the last terms in Ch,s describe the on-site repulsions, which 
soften the hard-core boson condition. It is noted that this effective theory is underpinned by a so-called bosonic RVB 
order parameter: A" = ^^^(^"^"^"^i-i-aCT^I-o-) 0, which describes a short-ranged spin liquid state in the spinion 
background in the underdoped regime of a doped Mott insulatoil^. 

The gauge fields, A'^ and A*, capture the essential physics of the phase string effeclPin a doped Mott insulator via 
the following lattice mutual Chern-Simons term in Eq. ([s]) 

Cues = ^ ef^"^ (A^ - 27r^/) (A^ - 2^^'') , (6) 

with {.vK"'''} the integer fields due to the compactness of the lattice field theorjEE^ which are important in restoring 
the spin SU(2) symmetry. Without j^"''*}, Eq. ([l]) is recovered the present theory reduces to the original mutual 
Chern-Simons gauge theory formulation of the phase string theory in Ref. |H1 
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B. Low-energy effective theory for superconducting phase 

Superconductivity can naturally arise in a saddle-point state of the above mutual Chern-Simons gauge theory, 
in which holons become Bose-condensed and spinous are gappecP^Sl. in the folUowing, we shall derive a general 
effective theory governing the low-energy fluctuations of the gauge fields A'^'^ around this saddle-point. This new 
effective theory will replace the conventional Ginzburg-Landau theory to describe the generalized rigidity of a non-BCS 
superconducting state. 

By minimizing the action S in Eq. ^ with regard to Aq and Aq, one can find the saddle-point solution of A'^'^, 
i.e., A'''* defined by 

e'^f^d^A'l = nS, = 0, (7) 



in the presence of a uniform distribution of the holon condensate with 6 as the holon concentration. One can always 

choose the temporal components = by shifting A'*^^ in Eqs. Q and (jsj). At the saddle point, one has .yf — 0. 
Equation ([t]) can be understood as that a fictitious uniform flux ttS is generated by the average distribution of the 
holon condensate via A^, while the vanishing net flux in A" is due to the fact that the spinous are short-range paired 

with (j2a '^^\<j^ia^ — (where (...) denotes the expectation value at the saddle-point). 

Around the saddle-point, the gauge field A^ may be rewritten as 

4 = ^'+<- (8) 

Then an effective low-energy Lagrangian of the "weak" fluctuations of and A^ around the saddle-point can be 
derived as follows 

£eff[^^«^ A^] = ||(^^ + + ^ {e'^'-'d^alf + ^e'-^A^a.a^ . (9) 

Here the first term on the right-hand-side (rhs) of Eq. ([£]) arises from Ch in Eq. Q after considering the holon 
condensation"^, with = tt^/ui and Qi — 92 — '^'^^Phth- Note that the gauge field A* is fixed in the unitary gauge 
by absorbing Goldstone boson field in the holon superfluid h — y/p^^e'^, and the superfluid density fluctuation 
Sph — Ph — Ph has been integrated out around Ph = S. 

The second term on the rhs of Eqs. ^ originates from the spinon term Cg in Eq. ([s]) after integrating out the 
spinon field b, with = ffi = 51 = (tli6 details are outlined in Appendix fAl) . Here, nis denotes the spin mass 
gap which is an important parameter as discussed in Appendix A below Eq. ( Af[). The parameters 7 > and k > 



as given in Eq. (A61. Note that in Refs. 16171 by assuming that the state is deep inside the superconducting phase, 
the limit of nis — > 00 has been taken such that gfi = there. Equation (joj) corresponds to the general case for a finite 
spin mass gap nis. 

The last term on the rhs of Eq. ^ is the mutual Chern-Simons term with the high-energy fluctuations of ./K*''' 
omitted. Note that the role played by in constructing a topological object like the magnetic fiux quantization 
has been discussed in Ref. [7] and the references therein. 

C. Meissner effect 

To see that the effective Lagrangian Ccs in Eq. ^ indeed describes a superconducting state, one may further 
integrate out A'^ and a'* to arrive at a low-energy action 



(10) 



where the summation is replaced by integration in the continuum limit, with the external magnetic field A*^-"- 
fixed in the Coulomb gauge (V • A*^^ = 0) . Here the transverse photon mass Mp is given by 



A/p - f (11) 



with the London penetration depth 



\p ^ J + KTUs^ (12) 
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in which An = ,„} - ■ 

Therefore, so long as rus 7^ 0, the system does exhibit the Meissner effect with a finite Ap in response to an external 
magnetic field. The physical electric current J*^ satisfies the London equation: 

- SA-^ "Ap ■ ^^^^ 



Equations (111 and (12) will recover the previous result!^ obtained in the "ideal superconductor" limit of — > 00, 
where the London penetration depth Ap approaches the "shortest" Aq. Namely, in the large spin-gap limit, the 
superfluidity of a holon condensate is well protected. 

However, in the opposite limit of — > 0, i.e., with the spin gap closing up , one finds Ap — > 00 in Eq. (12) such 
that the Meissner effect disappears. It implies that the superconducting phase coherence can be destroyeoby the 
gapless spin excitations at zero temperature, even in the presence of a bare holon condensate with a finite Aq. This is a 
peculiar manifestation that the neutral spin degrees of freedom can strongly influence the superconducting condensate. 
The quantum phase transition of superconductivity occurring here is quit different from the conventional transition 
of BCS-type as in the latter a charge (superconducting) gap must be closed up at the Fermi surface. Physically, with 
rris — > 0, the proliferation of neutral spinous will disorder the superconducting phase coherence via the 7r-vortices 
bound to them in the mutual Chern-Simons gauge theory. Such a superconducting phase transition can also occur at 
a finite temperature at T = "^s which is of non-BCS-type (see Ref. [T7)) . 

Finally it is noted that rn^ = with a spinon condensation is a necessary, but not a sufficient condition for the 
emergence of a true AFM ordering in the mutual Chern-Simons gauge theory. A so-called Bose insulating phase (BI) 
with nis = but Aq 7^ has been proposecP^ to appear as an intermediate phase between the superconducting and 
AFM phases. In other words, the present superconducting phase can be terminated at a critical point defined by 
rUs — >■ before a true AFM order can set in with further reducing doping concentration near the half-filling. 

III. TOPOLOGICAL CHARACTERIZATION OF SUPERCONDUCTIVITY 

The effective low-energy Lagrangian Cqs in Eq. ([9| has replaced the conventional Ginzburg-Landau equation to 
describe superconductivity in the mutual Chern-Simons theory. Now we examine the following question: What is the 
essential physical quantity that captures the nature of superconductivity, since the BCS order parameter (as governed 
by the conventional Ginzburg-Landau equation) no longer explicitly appears in the low-energy effective theory? 

To be sure, the off diagonal long order (ODLRO) of superconductivity is still described by the Cooper pairing order 
parameter in the phase string theory of doped Mott insulator^. Nevertheless, since the Cooper pairing amplitude can 
persist well beyond the superconducting phase, as determinecP by the bosonic RVB parameter A'' 7^ and the holon 
condensation (h) 7^ 0, the low-lying excitations and the phase transition of superconducting phase are essentially 
characterized by the topological phase fluctuations via the mutual Chern-Simons gauge fields, A'^ and a'', in the 
low-energy effective Lagrangian ([9|. 

Physically, A" and A^ are rel ated to the spin and hole currents in the mutual Chern-Simon gauge theory by the 
following equations of motiorPl^ 

j-=leO"^iJ|, j^ = ie°"^4\ (14) 

where the internal "electric" fields in the imaginary time representation E^'^ = — ie^^^d^Af,''' and the spin and charge 

currents are defined as j^^^ = f^h/l ■ Hence the low-lying fluctuations of A" and a^ in Eq. ^ actually describe the 
current fluctuations of the spin and charge degrees of freedom, while the amplitude fluctuations of the matter fields 
remain gapped in the superconducting phase. 

A. Wilson loops 

In order to characterize the gauge fiuctuations of A^^'*, a pair of Wilson loops have been introduced^ 

W''-''[C]ee (w'-''[C]'j , (15) 

where 

{...)=Z-^ J D[A',A'']e-^''^''"... 
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and 

^^■"[C] ^ exp [iY,A^\x)J^{x)l (16) 

X 

in which the (2+l)-diniensional unit current, for a test holon or a test spinon, Jq{x) = +1 (—1) for the Unk x ^ x + jl 
(x + /x — >■ x), and is zero otherwise on a close loop C. 

Physically, the Wilson loop VF*'[C]or 14^'' [C] probes the interaction of a pair of test holons or spinons at a distance 
R, defined by V''{R) or V'{R), vifP 

V'^-{R)^- lun ^InW^^'n (17) 

T— >-oo 1 

if C is taken as a spacetime rectangle with length T (R) in the imaginary time (spatial) direction. In Refs. 16171 
in the limit of large spin gap {nis — >■ oo), one finds V {R) ^ and y*(i?) ~ Ini?. Namely the holons are 
deconfined while the spinons are logarithmically confined. It is consistent with the present saddle-point state of holon 
condensation discussed in Sec. II. 

In the presence of a finite spin gap to^, we find that the above conclusion remains true, with the the confining 
potential, between a test spinon and test anti-spinon, modified to 

IT (R 



^^(^) = ^ln(-) (18) 



where Aq is replaced by Ap defined in Eq. (12 1 [R* is an ultraviolet cutoff of the distance). The result indicates 
that the spinons will experience a confinement-deconfinement transition at Ap — )■ oo or rus — > 0, coinciding with the 
superconducting phase transition critical point at T = 0. On the other hand, V^{R) for holons remains deconfined. 



B. Topological characterization 

The Wilson loops, iy''*[C], can provide important information on the confinement/deconfinement of spinons and 
holons in the superconducting phase as discussed above. However, to understand why W^*''^[C] behave as such, the 
topological nature of the mutual Chern-Simons effective theory in Eq. ([9| has to be further revealed. 

For this purpose we define 

e'^ = K,-^ (w'lCi] W^\C2]) , (19) 

with the normalization factor 

JC = W^[Ci]W'[C2\. (20) 

Let us first consider the case of a pure mutual Chern-Simons term Q without coupling to the matter fields [i.e., 
omitting the first two terms in Eq. ([9|], a straightforward calculation shows 

I = 7re (21) 

where G = 0, ±1, ±2, • • • , denotes the Gaussian linking number, which is the winding number of the close loops, Ci 
and C2, as shown in Fig. [l] (a) and (b). To obtain Eq. (21 ), we write down the explicit expression of X in a pure 



mutual Chern-Simons gauge theory. 



where, the unit spacetime currents J''^ and J''^ keep track of spacetime loops, Ci and C2, respectively. Integrating over 
A** leads to a flux constraint on a'' configuration, namely, e^'^^dvofl = tt . It means that the gauge field a'' forms 
a spacetime closed 7r-fluxtube along Ci. / d'^xa^J^^ = af^dx^ denotes the spacetime flux piercing C2. Noting that 

there are Q times for the fluxtube Ci piercing C2, one arrives at e*-^ = e''^®, with /C = 1. Ignoring the 27r-periodicity, 
we thus obtain Eq. (21). In Witten's seminal papei^, a similar Gaussian linking number structure was found in the 



Chern-Simons gauge theory, which is related to the Jones polynomial in the Knot theory. 
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(a) G=l (b) G=2 

FIG. 1: (Color online). Schematic picture of two loops Ci and C2 in (2+l)-dimensional Euclidean spacetime. Ci is an arbitrary 
loop but restricted inside spatial plane (i.e., the xy plane) at some imaginary time. C2 is an arbitrary loop with a definite 
Gaussian linking number G Z (0 = 1 in (a); = 2 in (b)) that characterizes the number of times that C2 winds around Ci. 
The ball (in blue) and the arrow (in red) stand for the test holon and test spinon, respectively. The arrows on the two loops 
indicate the direction of test spinon current and test holon current. 



In Fig. [2] a gedankenexperiment is shown to illustrate the case of winding number = 1 in the (2+l)-diniensional 
Euclidean spacetime, with a pair of test spinon and antispinon creating and annihilating on a close loop C2. If the 
background spin degrees of freedom are gapped with — > 00, the test spinous cannot induce spin excitations from 
the ground state and thus the net contribution to e*-^ comes from the tt flux phase shift (with 9 = 1) mentioned above 
due to the pure mutual Chern-Simons action. 

On the other hand, if the spin gap becomes finite, then the spin excitations can be induced from the spinon 
background to renomalize I. Now we consider the full £off given in Eq. ^ . By choosing the close loop Ci of the test 
holon within the spatial xy-plane, one finally arrives at the following compact result (for the details, see Appendix 
B): 

I^irQ {Xo/Xpf . (22) 

We see that I remains to be real and positive as a topological quantity, which is invariant under a continuous 
deformation of loops Ci and C2 for a given Gaussian linking number Q, as illustrated in Fig. [T] But I is renormalized 
by a factor (Ag/Ap) which is solely determined by the spin gap m^. In particular, in the limit of — > 0, one finds 
1^0. 



If one rewrites Eq. ( 22 ) by 

Ap = Ao • (23) 

where Ap becomes a function of topological order parameter X for a given Gaussian linking number = 1, then the 
limit Ap — > 00, i.e., the disappearance of superconductivity, will correspond to I = 0. In other words, the topological 
quantity X can be regarded as the superconducting order parameter in the effective theory of Eq. (|9|, and the quantum 
phase transition at T = occurs when I gets fully screened by the gapless spin excitations at rus = 0. 

As illustrated in Fig. [l] Z 7^ indicates that a test spinon can be always "felt" by a vortex-like current loop of a test 
holon if their mutual winding number Q ^ Q. Physically, since the holons are in the Bose condensed state, a spinon 
excitation will then induce a supercurrent vortex response from the condensed holons throughout the sample to form a 
so-called spinon-vortexi^^. This is the origin of spinon confinement discussed above, since each spinon-vortex composite 
is logarithmically divergent in energy which will result in vortex-antivortex confinement. In a non-superconducting 
phase where the spinous become deconfined due to the vortex-antivortex disassociation and proliferation, one finds 
X = in the gedankenexperiment shown in Fig. [2j Namely, 2 = describes a phase disordered superconductor. 



C. Experimental prediction 



We propose an experimental setup to probe the topological characterization I of the superconducting phase in Fig. 
Isj Note that 1^0 reflects the mutual statistical interaction between spin and charge degrees of freedom (cf. Fig. II]). 
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(a) x=xi 


(b) x=X2 








■ * 


(c) X=X3 


(d) x=X4 



FIG. 2: (Color online). Gedankenexperiment (Gaussian linking number p = 1 for simplicity), (a-d) are four snapshots of spatial 
plane subsequently at imaginary time t\ < T2 < tz < t^. (a) A pair of test spinons with opposite spin indices are created at 
a spatial site "P" at t\. (b) Two spinons arrive at new spatial positions at ti via the dashed trajectories, (c) At rs, a holon 
current circuit Ci is located in plane, indicating that an Aharonov-Bohm phase around Ci is picked up by the test holon. (d) 
At r4, the pair of test spinons meet each other again at some position and then are annihilated into vacuum simultaneously. 

By using the idea illustrated in the gedankenexperiment of Fig. [2] one may inject free "test" spinons by heating up 
the middle semi-island part of the sample in Fig. [3] to a temperature slightly higher than T^, such that free spinons 
can be excited and flow towards the lower T end connecting to the outside superconducting ring. With T 7^ inside 
the ring, one expects that such an injection of free spinons can generate supercurrents in the superconducting ring in 
a form of significant noise spikes due to the unpolarized spinon-vortex effect. 




FIG. 3: (Color online) A proposed experiment to detect mutual statistics in a mutual Chern-Simons superconductor. The 
temperature in the middle semi- island is raised slightly above the superconducting critical temperature T^, while the outside 
part remains below Tc and spinons are confined into vortex-antivortex pairs. The confinement is indicated by the black wavy 
line between two spinons. The small directional discs (red arrows) represent vortices (spinons). In such a mutual Chern-Simons 
superconductor at T > Tc, spinons become deconfined to form free spinon- vortices^. Then, the fiow of such free spinons 
from the semi-island towards the cooler end connecting to the outside superconducting ring is expected to generate random 
supercurrents in the latter as indicated by J*^ (blue dashed line). 

Instead of raising the temperature, one may also apply an in-plance magnetic field to reduce the spinon mass 
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gap nis inside the semi-island part by = rus — ^gLf^sB'^, where ql is the Lande factor and is the Bohr 
magneton. Then a sufficiently strong B'^ may drive the semi-island part into a spinon-vortex liquid state with X — 
via the Zeeman effectJ^^. At the same time, the gradient of can also naturally cause the free spinous flow into the 
superconducting ring to induce the similar supercurrent response shown in Fig. [3j 



D. Emergence of fermionic excitations: A proof-of-principle derivation 

The mutual Chern-Simons theory given in Eq. (|3| is an all-boson description of doped Mott insulators. In the 
above subsections we have studied the low-lying collective fluctuations of this theory based on the effective Lagrangian 
([9]) . An important question that arises naturally is if and how a low-energy fermionic excitation can emerge in such 
an all-boson theory, which is independent of the above collective excitations. 

To give a proof of principle, we first note that the mutual Chern-Simons term in Eq. (|9| can be replaced by an 
auxiliary Lagrangian in real-time representation: 

/:aux = (^/*^*/ - A*/*/) , (24) 

where is an auxiliary 4-component Dirac spinor in (2+l)-dimcnsional spacetime with "flavor" indices / = 1,2. 
'P = ^^'Dfj_. Here 7-matrices are defined in 4-dimensional reducible representation of 2-dimensional Clifford algebra, 
i.e., 7*^ = diag{a^ , —u^) , 7" = diag{a'^, where a = 1, 2 and cr^'^''^ are Pauli matrices. The covariant derivative 

is defined as = 9^-1- iAJ^ + ir^a'^ The 4x4 matrix T3 = diag{I, —I), where I is the 2x2 unit matrix. The 



second term in Eq. (24 1 is mass term with A > 0. 

£aux preserves a U x U{1)v gauge symmetry where subscripts A and V denote axial vector type and ordinary 
vector type U{1) gauge invariance, respectively. It originates from the fact that is an ordinary vector field with 
gauge group generator "1" while is an axial vector field with gauge group generator t^. The parity and time- 
reversal symmetries are conserved even in the presence of finite mass gap which may be expressed as: A^I'j'I'j — 
A*y7o*/ = A*J,^cr^vl/^^ _ A*J,_(T^^'/_, where, (*/-) is upper (lower) 2-component spinor of Here, ± 
denotes the chirality^". In the case of parity (P is parity operator). 



It is easy to check that the mass term is invariant under such parity transformation in (2+l)-dimensional spacetime. 
Integrating over ^ f leads to a pure mutual Chern-Simons term by ignoring Maxwell terms (oc momentum squared) in 
the one-loop quantum correction (Fig. |4]) in the long wavelength limit \p\ A. The key s teps o f derivation are given 
as follows. As shown in Fig. [4] the vacuum polarization for and a'* at one-loop leveP^ZSl ig a 3-dimensional 






FIG. 4: One-loop quantum correction. Zigzag, wavy, and double lines represent ^4", a'', and 'I', correspondingly. The first two 
diagrams correspond to Maxwell terms, while the last one gives the non-vanishing mutual Chern-Simons vertex. 

Euclidean momentum): 



n(p) - -^^^ 



■ p2^4A=^ . f p 
2A -I arcsm 



47rp2 ' p V(P^ + 4A2)i/2 

where the total flavor number A/ = 2. The expansion around small p is given by 

AO' AO'P^ 3A/>p4 



H(p) 



67rA 607rA3 11207rA5 



9 



if" f 

The propagators of gauge fields can be constructed in the Landau gauge and the Maxwell terms read out: — 12^ a — 

j^2^^ with f^j^ = 9^ A* — and fj^^ = S^ajj: — d^a'^^^. Meanwhile, the nonvanishing mutual Chern-Simons term for 

A* and a'' in Fig. [4] is also radiatively generated without breaking parity and time-reversal symmetries. Note that, 
upper(lower)-component spinor has mass A (—A): 



(27r)3(fc2 + A2)((p + A:)2 + A2) An ^"^"^ Jo ^A^+£{l~£)p^ 

where, the superscript a =t denoting upper and lower 2-component spinor contributions respectively, and 

Ifi = jjc^i cr^, tj^} with identity: Tri^^^vlx) = 2e^yA- In the sufficiently low momentum transfer limit {p — )• 0), it is 
expanded as power series of p: 



In Eq. (261, the sign factor "cr" renders vanishing of Abelian Chern-Simons terms (~ e'^'^'^'^A^'^d^A^ and ^ e'^'^^a'^d^a'^). 
But due to the structure of matrix ra, the crossing term, i.e the mutual Chern-Simons term survives as expected. 
Mathematically, the above quantum corrections are organized as below: 

i ^ (A^ + aa';,) ^e^-'d, {Ai + aa^) ^ le^^^'A^d^^ . (27) 

where the Maxwell terms are ignored in the long wavelength limit. 

It is important to note that the fermionic auxiliary field 5*/ is not gauge-neutral as it still couples to and in 



Eq. (24). A true fermion excitation of quasiparticle should be gauge-neutral without knowing A* and a^. It is only 
possible if there exists a residual attractive potential such that a holon and a spinon excited from the condensates can 
form a bound state properly with the auxiliary field f . The resulting object ^ will be both fermionic and coherent 
because all statistical gauge charges carried by ^/ can be exactly canceled by that of the spinon and holon partners. 
In Fig. [5] the spacetime loop of a composite-particle ^ is illustrated as a gauge neutral and thus physically observable 
elementary excitation at low energies. 

In particular, it carries the right quantum number as a Bogoliubov quasiparticle. As noted before, to produce the 
right prefactor j/tt of the pure mutual Chern-Simons term, ^I*^ has totally 2x2 = 4 branches of positive energy 
solutions, which are consistent with the existence of totally four quasiparticle nodal regions in the Brillouin zone, 
although the detailed energy spectrum depends on the residual attractive potential which goes beyond the present 
long-wavelength effective description. On the other hand, the existence of a nodal fermionic Bogoliubov quasiparticle 
has been recently demonstrated based on a ground state wavefunction approach to the t-J model in the phase string 
formulation, in which the quasiparticle as a similar bound state has been explicitly shown^. Therefore, the present 
field-theory formulation provides a proof- of-principle for a fermionic quasiparticle mode emerging in the mutual 
Chern-Simons theory. 




FIG. 5: (Color online). 'I' is composed by a spinon (red arrow), a holon (blue ball) and a Dirac fermion (small grey ball ^ f). 
Its spacetime loop is composed by three loops rendering gauge charge neutrality. 
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IV. CONCLUSION 



In this work, we have obtained an effective theory of the low-lying fluctuations in the superconducting state of the 
mutual Chern-Simons theory. In contrast with a conventional Ginzburg-Landau action for a BCS superconductor, 
the present effective description in Eq. ^ deals with peculiar phase fluctuations in the doped Mott insulator via a 
pair of mutual Chern-Simons gauge fields, which are related to the fluctuations of charge and spin currents of the 
condensate. 

The superconducting phase coherence is found to be properly characterized by a topological order parameter which 
measures the mutual entanglement of the charge and spin currents in the low-energy effective theory. In particular, a 
zero-temperature quantum phase transition is shown to occur, due to the vanishing of this topological order parameter, 
when the spin gap is reduced to zero at low doping. 

We have discussed a possible experimental probe into the long-range entanglement between the supercurrents and 
the neutral spin excitations as a unique signature of this non-BCS superconducting state. We have also given a 
proof-of-principle how a fermionic quasiparticle can emerge in the mutual Chern-Simons theory. Together with the 
previously identified flux quantization^ * at hc/2e, the non-BCS superconductivity in the mutual Chern-Simons theory 
is well established. 
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Appendix A: Path-integration over spinons 



First of all, we note that the spinon Lagrangian Lg m Eq. ^ may be further siraplified to a CP(1) form m the 
long-wavelength, low-energy regime as followP 



\{df,-iaA^^)z„\ +ml\z„\'^ 



(Al) 



Here the temporal component is rescaled by the spin wave velocity Cg = \/ Jsi^s + 4Js) with the coupling constant 
g — Cs/ (2Js), and is treated as a dynamically generated spi non mass enforcing l-Zo-P — 1- 

The problem of studying the fluctuation of ^4^, or a^, in Eq. ( Al ) is then transformed into calculating the vacuum 
polarization tensor in (2-|-l)-dimensional scalar quantum electrodynamics in the presence of a uniform background 
magnetic field of A'* defined in Eq. Q. The technical challenge arises from the fact that because of the momentum 
is no longer a good quantum number. Nevertheless the issue is resolved by using the so-called "string-inspired" 
formalisirPS. Following the procedures and conventions in Ref. |24[ we obtain the polarization tensor: 



n^^(fc) 2 X 



(47r)^ 



~T 



2 e 



-™'^det-5 



"sin(Z)' 




Z 


I 



(A2) 



where the prefactor 2 comes from the summation over the spin index, T denotes the Schwinger proper time, and 
spacetime dimension D — 3. Here Z is a 3 x 3 matrix, defined by Z = nFT with 



(A3) 




The determinant in Eq. ( A2 ) can be simply factorized into 

, _i fsinfZ) 
det . 

and the transversality of Lorentz tensor lui, is given by 



■ST 



sinh(7r5 T) ' 



_ sinh(7r^T(l-2£)) _ 



(A4) 



(A5) 



By further restoring the spin wave velocity c^, we arrive at the "Maxwell" term on the rhs of Eq. ^ which describes 



the fluctuation a!^, with the coefficients given by 



7 = 



12c,, 



12 



(A6) 



in which the dimensionless coefficient i? depends on the dimensionless quantity 



as follows 



/o°° ^^e-"^ /q d£ (1 - 2^) • TTc, S Tmj^ sinh ((1 - 2£)7rc, 5 Tmj'^) ■ (sinh(7rc, 5 Tmj'^)) " 



Jo^^e-^j^deii-2iy 

Note that ■& — 1 in the limit of the absence of "magnetic flux" in A^. 



(A7) 



Appendix B: Derivation details of topological order parameter X 

To obtain I, we first note that the Lagrangian with holon and spinon currents circulating on loops Ci i (m* = 



12 



is quadratic in the gauge fields thus allows rigorously integrating out A'^ and in the generating functional 

(W' [Ci]W'' [C2]) = J i^[AV]e-J''^'^^^cP„c.] ^Const. x J D[a''] e^v J {e^,xd,a\f ~ la^^J^ 

(B2) 

in which we have defined 



^ \ J 

and absorbed terms unrelated to mutual correlation between Ci and C2 into Const, which will be canceled by the 
denominator in the definition of I. 

Adopting the "deformed" Lorentz gauge condition, 

T^SoaS + ^dia'l + ^32^ - (B3) 

(_ro Lrl (_r2 

and integrating out in the Feynman gaug^^, we get 

iw' [Ci]M/" [C2]) = Const. X expj-i^-^^ j Sx J d^y ^ G,, J^(a;) J^(y)/(x, y)| (B4) 

in which 

-1 



Expanding the integrand in Eq. |B4| and extracting mutual correlation terms between Ci and C2 only, and noting 
that / d^xJ^^{x) — §^ dx^, i — 1,2, we get 

+ dyi(^2<j> dx2dlf{x,y) - Go j> dxody{x,y)^ (B5) 

+ ^2/2(^00^ dxoaf/(x,?/) - Gi^ ^ 

Applying Stokes' formula to transform the loop integral over Ci into a surface integral over a disk Di such that 
dDi = Ci, and noting that dx^d1^h{x) = 0, we get 

+ f dt„AdvJ dx,Gj ^(-u,kl-'^k})(j2^] •r"-"-" (B6) 



(2^)3 V ' ^ V V A 

+ ^ dy2Adyo£ dx^Gi j (""y ^0 " "1^2) 



Given that Gi = G2, Eq. B6 is greatly reduced if Ci lies within a spatial plane, 

-1 



1 TT^ 2 



/ dy, A d,2 / dxo^Go / ^ (-fc? - .2^ - ^ki) E 



1 g-ii-^-y)-k 



2GaGiG2^JD, Jc, 5 "7 (2^)n ^0 /Va^\ 

1 TT^ 2 m* /■ /■ 

777^7^ ^GoGi / dyi A fi2/2 $ fixo-J^^) [x - y) 

1 G0G1G2 TT d J^,^ 

(Ao/Ap)Ve 
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in which Q is the Gaussian Unking number of Ci and C2 loops. To obtain the last line of above identity, the definition 
of Xp and Aq in Eq. ( 12 ) is applied. 



